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Abstract 

Spatially homogeneous and isotropic cosmological models, with a perfect fluid matter 
source and non-vanishing cosmological constant, are studied. The equations governing 
linear perturbations of the space-time and the variation of energy density are given. The 
complete solution of the problem is obtained for C°° perturbations, using a comoving time. 
The Sachs- Wolfe fluctuations of the temperature of the cosmic background radiation are 
obtained for the relatively growing density perturbations. It is found that the observable 
celestial microwave fluctuation pattern underwent a reversal approximately two billion 
years ago. What is observed today is a negative image of the last scattering surface with 
an attenuation of the fluctuations, due to the presence of the cosmological constant. 
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1 Introduction 



Recent observations by the High-z Supernova Search Team ( Tonrv et al. . 2003f ) corroborate the 
data on an accelerating expansion of the Universe. Three independent lin es of evidence [those 
from the Wilkinson Microwave Anisotropy Probe (WMAP) measurements ( Spergel et all 12003) 
of the cosmic m i crowav e background radiation (CMBR), the Sloan Digital Sky Survey (SDSS) 



(Te gmark et all l2004al ) and the observations of type-la supernova spectra] converge on the 
value Q\ = 0.70 ± 0.04 of dark energy. All data are consistent, within a 20% error bar, with 
a time-independent dark energy distribution, as is the case with the cosmological constant A 
and flat space (k = 0). These developments attracted attention to cosmological models in the 
presence of a A term in Einstein's gravitational equations. Many excellent revi e w papers are 
availa ble (|CarrollL 1200 it IPeebles fc Ifatral 120021 : ICarroll. Press, fc Turner! . Il992t IWhite et al. 



19941 1 on the implications of a cosmological constant in models of the Universe. While the 
cosmological constant has little effect on the large-scale structure and dyna mics of the Universe , 



the microwave fluctuations will be significantly different in A ^ models ((White et allll994h . 



The scalar perturbations of A ^ models with multicomponent fluids (i.e., cold or hot dar k 



matter, photon s and massless neutrino) have been com p uted by Bond fc Efstathio"ul (Il987 ). 
Fukugita et al. (|l990i) . iHoltzmanl (|l989lb IVittorio et all (|l99lh . iHu fc Sugivamal (|l995h and 
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StomDor (1994) in the linear approximation. These numeri cal results have b een c ompared by 



1994 ) and found to be in reasonable agreement. Hu fc Sugivama ( 1995j ) claim that 



for A ^ 0, one must use a numerical approach. In our work, however, we present an analytic 
treatment of the collisionless dust which represents the late-time evolution. 

A complete solution describing the effect of the cosmological constant on perturbation 
dynamics is lacking in the literature, this including the rotational and gravitational wave 
perturbation s. The temperature fluctuations of the CMBR have been computed in numeri- 
cal schemes dBond fc Efstathiou L 1987 L Fukugita et all Il990l: Holtzman , 19891: Vittorio et al 



ll99lUStomD'orr Il994j ) and lMultamaki fe Elgarovl (|2003| ) investigate the integrated contributions 
Known as the Sachs- Wolfe effect, this is the variation in the redshift of a photon travelling freely 
in the universe rippled with perturbations. Initial work, both on the zero-pressu re model of the 
curren t state of the Universe, and on the radiation-dominated era, was done by Sachs fc Wolf el 
f 19671 ). The pressure -free universe in the pre sence of a cosmological constant has been consid- 
ered by Heath (1977) and Lahav et al.l (|l99l| ). They used the re dshift paramet e r to in vestigate 
the behavio ur of the density cont rast. A more recent work by Vale fc Lemosl ( 2001 ) uses the 
formalism of Padmanabhanl ( 1993) to compute linear perturbations of a dust-filled Friedmann- 
Robert son- Walker (FRW) universe in the presence of a cosmological constant. They find that 
the presence of A inhibits the growth of the fluctuations in time. [A different approach to the 
perturbation proble m of a different (de Sitter) model, using the gauge-independent formulation, 
has been pursued by iBarrow (2003). for the investigation of stability under perturbations.] 

Neither of these earlier works reaches a complete perturbative picture. Although the Sachs- 
Wolfe effect is computed by numerical methods, there is no analytic treatment so far available 
for A ^ 0. Our aim here, therefore, is to consider, in an analytic framework, the Sachs- Wolfe 
effect and its contributions to the fluctuations of the CMBR i n the presence of a cosmological 
constant. We follow the treatment of the perturbations due to White! f 1973! ) . thereby relaxing 
the momentum conditions of the origi nal Sachs- Wol fe work. Both works yield all C°° perturba- 
tions although this was only shown by White! ( 19731 ). Correspondingly, in the present paper we 
also obtain all C°° perturbations. In Sec. 2, we enlist the perturbed field equations using the 
conformal time rj. While retaining the form of the metric functions, we change to the comoving 
time t in the field equations so obtained. The reason for this is that the unperturbed radius a 
has a simple analytic dependence on t. 

The general solution of the linear field equations in a dust-filled universe is obtained in Sec. 
3. We find that all first-order fields, except the wave solutions, can be expressed in terms of 
only two complex incomplete elliptic integrals E and F. This is achieved by using the mirror 
symmetries of the elliptic integrals. 

In Sec. 4, we obtain the Sachs- Wolfe effect and its integral contribution which vanishes in the 
A — > limit. An unexpected feature of the transfer function is that it changes sign and reverses 
the pattern of the temperature fluctuation on the celestial sphere at late times. In Sec. 5, we 
investigate the physical interpretation of our solution. Taking the current experimental values 
of the cosmological parameters, we obtain the comoving time elapsed between the emission 
at the surface of last scattering and reception of the photon. We find that the Sachs- Wolfe 
contribution to the CMB fluctuations in the presence of the cosmological constant is damped 
by a factor of 3 today, because of the proximity, on the cosmological time scale, of the moment 
of reversal of the microwave temperature fluctuations and the reception. This feature of the 
SW effect has remained concealed in earlier works, because of the numerical methods involved. 
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2 Linear perturbations 



In this section we present the equations of a perturbed spatially flat (k = 0) FRW cosmology 
in the linear approximation. The metric is that of a perturbed FRW model in the conformal 
form 

g a b = a 2 (r])(r] ab + h ab ) , (1) 

where r\ a b = diag(l, —1, —1, —1), h a b is the metric perturbation and the scale function a = a{rf) 
is determined by the field equations. Roman indices run from to 3 and Greek indices from 1 
to 3. Except when we compute numerical values, we employ units such that the speed of light 
is c = 1 and the gravitational constant G = l/8n. We use the Minkowski metric T) a b and its 
inverse to raise and lower indices of h a b and of other small quantities. 
The matter source is assumed to be a perfect fluid, 

T ab = {p + p)u a u b - pg ab (2) 

with the four velocity u a of the fluid normalized by u l Ui = 1. The cosmological constant A 
of the Einstein equations is introduced in the energy-momentum tensor. The effective energy 
density p and pressure p are related to the corresponding fluid quantities as 

p=p+A, p—p—A. (3) 

In the perturbed space-time we choose comoving coordinates ( Sachs fc Wolfd . 1967 ). Using this 
gauge, the four velocity of the fluid and the metric perturbation h Q0 are 

u a = ^ , h 00 = . (4) 
a 

The coordinate transformations x a — > x a — £ a , which preserve the coordinate condition, have 
the following properties 

Cu a r^a Q ^ e = K^) ; C = c « {x P h (5) 

a 

where the arbitrary functions b and c a are independent of the conformal time r). 
The conservation of energy-momentum 

T ab = (p + dp + p + 5p)u a u h -(p + 5p- A)g ab (6) 

implies the unperturbed relations 

3(p + p)a' + ap' = , p, a = (7) 

and the perturbative equations 

5p' + 3^(8p + 5p) + (p + p)^ = 0, (8) 
(h a0 ' + ^h a0 ^j (p + p) + h a0 p' -5p a = (9) 

where h = r] rs h rs . The prime denotes derivative with respect to the conformal time r\. 
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Einstein equations G a b = T ab to leading order give 



■A 



P + A 



a" a' 2 
> 

'a 3 a 4 



-p + A . 



(10) 



The perturbed field equations for a perfect fluid have been obtained bv lwhitel (jl Q73h . They 
remain valid in the presence of the cosmological constant provided the pressure and density are 
replaced by their tilded values according to ©• The scale factor a is now a solution of Eqs. 
()10|) . Combining the 00 component and the trace of the spatial part of the Einstein equations 
we have 



h " = -2-ti + 2h°\' + 4-h°\ 
a ' M a 

+\S^ U + l -V 2 h - 3a 2 5p 



S^u + ~V 2 h + 2- (2h\ - ti) + 2a 2 5p = . 



;i2) 



In these equations S a p = h a/3 — rj a/3 h/3 denote the trace- free part of the perturbations and 
V 2 / = —rj^f^u, where V 2 = A is the standard Laplacian, for an arbitrary smooth function /. 
The remaining components of the linearized Einstein equations are 



V 2 h 0a - -h a ' + S°v ' 
3 ' ^ 



'a" a! T 



,0a 



(13) 



sy+2^sy-v 2 s a p 



+ h a °' + h, " 



a 



-V 2 /i + 2/i°V + 4-/*°\ 
3 ,A4 a * 



(14) 



3 Integration of the field equations 

In this section we solve the field equations in the presence of pressureless matter, p = 5p = 
for perturbations in the class of C°° functions. We then repeatedly use the following 
Lemma (Brelot) 

If g is any C°° function on E 3 , then there exists a C°° function f on E 3 such that V 2 / = g 



(A modern proof was provided by Friedman ( 19631 )). We introduce th e comoving time coor- 
dinate t by the relation adrj = dt. The solution of Eqs. (jlO|) for matter ( Stephani et al. . 2003) 



is 



p = , pa 3 = C M , a = a sinh 2/3 (Ct + C 
A 



(15) 



a 



C 



3A /Ax 1 / 2 

— ,H = coth( V3At ) , 



where Cq = Cto and Hq is the Hubble constant. With a new definition of t we set Co = 0, and 
the big bang occurs at t — 0. 
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The value of the deceleration parameter ( Misner . 19731 ) 



a /a 
a \a 



(16) 



decreases after the big bang, and changes sign (acceleration) at the moment of turnover t = tx- 
A dot denotes the derivative with respect to the comoving time coordinate t. The deceleration 
curve is displayed in Fig. 1. In the limit of a vanishing cosmological constant, the turnover 
time tx moves to infinity, and the radius has the time dependence a = dot 2 / 3 . 




Figure 1: The deceleration parameter q as a function of the dimensionless comoving time Ct. The 
present age of the universe is Ctn = 1.21 and the acceleration commences at comoving time Ctx = 0.66 
(Cf. Sec. EJ. 

The functional change of h a t> induced by the transformation (jSj) is expressed by the Lie 
derivative of the metric g a b with respect to £ a , 

^oo -> h QQ , h 0a -> h 0a + — , h -> h + 2c M + 66- (17) 

a ' M a 

ha/3 h a/ 3 + c a ^j3 + cp jCt + 2b-r] a!3 . 

(X 

Eq. (JHJ) i s integrated f or h 0a . With a coordinate transformation and using Brelot's Lemma, 
one can set (|Whitel . Il973h 

a 

where C a = C a (ar) is a spatial function with 

= 0. (19) 

As a consequence /i ^ = holds. Gauge transformations which preserve the above form of h 0a 
satisfy V 2 6 = 0. 

Integrating Eq. (jSJ) with respect to the conformal time we get 

h=-2— + H, (20) 
P 

where H = H(x a ) is an integration function. 
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From the perturbed field equations (jlljl and (J12)) we have 



h" = -2-h' + \s^ u + \v 2 h 
a 2 ,p 3 



-V 2 h 



From the above expressions 



h" + + a 2 5p = . 



(21) 
(22) 

(23) 



Inserting here h from Eq. (j2()j) . changing to the comoving time t and using the unperturbed 
solution (fTKj) . the equation for the density perturbation is 



a 2 C 



dp + 8-5 p + 3l-+4^--^)(5p = 0. 



a 2 Qa 3 



(24) 



Substituting the field equations (fTUj) we have 



+ 5 \8p = 



Changing to the variable 



results in the homogeneous equation 



z = cosh(2Ct) 



3(, 2 -l) (,-1) 



d 2 5p 
d5p 



(25) 
(26) 

(27) 



+ fllz 2 - 3z - 8) — - + (52 + l)5p = 
v 7 dz 



This is a special case of the Heun equation Application of the Kovacic 

algorithm ( Kovacic . Il98fil ) yields the solution 



Sp 



Vz + 1 

Z - l) 3 / 2 



K -k r dz ' 

1 2 A (z'-iywz' + i 



)3/2 



(28) 



where Xi = Ki(x a ) and = K2(x a ) are space functions. This contains the elliptic integral 

dz' 3 r x f x' 



h (^'-l) 1 /6(^ + i)3/2 2 2 / 3 y V^' 3 + l, 
where we have introduced the new variable x by 

z = 2x 3 + 1. 

Hence, using Eq. (J2Uj) 

h = H - coth(Ct) (2^ - K 2 J) . 



3/2 

I dx' 



(29) 



(30) 



(31) 



Taking the divergence of (fT3*j) . we obtain a total time derivative. By Brelot's Lemma, the 
integral of this equation defines V 2 B such that 



V 2 B = l -S^ v + l -V 2 h 



(32) 



where B = B(x a ) is a space function. From (jllj) we have 



h" = -2-ti + V 2 B . 
a 

Substituting here the form (|31j) of h, the terms proportional to K\ cancel and we get 

3V 2 B 



2^a 2 C 2 ' 



(33) 



(34) 



Since K\ is a C°° function, it may be represented ( Friedman! 1963) by a new spatial function 
A{x a ) as follows, 

Kt = V 2 A . (35) 

The time dependence of the density contrast 5p/p is displayed in Fig. 2. The contributions 
of the terms 

5 A = coth(Ct) , 5 B = coth(Ct)/ 

proportional to the spatial functions A and B, respectively, are relatively decreasing and grow- 
ing. They tend to constant values in the asymptotic future, 



lim 5 A = 1, lim 5 B = 2- 2/3 J Beto(|, |) = 1.086517930 
where Beta is Euler's integral of the first kind. 



2 

1.75 
1.5 
1.25 
8p/p 1 
0.75 
0.5 
0.25 




Ct 



Figure 2: The relatively decreasing and growing contributions to the density contrast. 
These terms have the series expansion in powers of r = Ct 

5A = - + ^-lr 3 + 0(r 4 ), 
r 3 45 

5 B = 2 1 / 3 -r 2 / 3 + 0(r 8 / 3 ). 



(36) 



In the limit of a vanishing cosmo logical constant, these are the density perturbations of Sachs fc Wolfei 

(23). 

If we rewrite Eq. (J13J) by introducing the new variable 



w = sinh(Ct) 
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we have 



d ( 2 \ V 2 V 2 C Q 

_ [ Q a L* h a ) j 



8C V 2 C a 



. 



3 w 2 y/l + W 2 
The integration of (|3*7jl and use of Eq. (|3*T|) yields 



5^ + -coth(Ct) (K.-K,!) ^ 



8CV 2 C a VI + w 2 



w 



v 2 v 2 c° 



+ 



dw' 



(37) 



(38) 



Here J a = J a (x^) are integration functions. From (fTTjl . the trace- free part has the trans- 
formation properties 



(39) 



The generator c a of the gauge transformations can be chosen (|Whitei Il973h such that the 
integration function in (f3*T]) is H = 3B and in (|3*%|) we set J a = 0. The normal form of the 
elliptic integral 

r w dm 1 

J= I — ^ — (40) 



w '4/3(l + 



W 



mV 2 



will be given in the Appendix. The graph of J is displayed in Fig. 3. The following relation 
holds between the time variables: 



x = — = sin 
a 



h 2 / 3 (Ct) 



2 ) 



1 



x 6 = sinh 2 (Ct) = -(cosh(2Ct) - 1 



z - 1 



2 V 2 
where the ranges are < t < oo, < x < oo, < if < oo, l<z<oo. 



(41) 
(42) 



-2.5 
-5 
-7.5 
J -10 
-12.5 
-15 
-17.5 
-20 



Ct 
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1.5 



Figure 3: Time dependence of J 



S 



Equation (|14|) for the trace-free part of the perturbations is 

a p + SaaS^ - W 2 S a p = h a p -AB a p 

5* ' 
3 



a 2 S a p + ZaaS a p - V 2 S a p = h a p - AB « (43) 
-2CV 2 {C% + C p a ) coth(Ct) + -^{V 2 h - 4V 2 S) 



V 2 V 2 (C^ + C A «) 
a 2 C 

with h as given in Eq. (|3*T|) . 

A particular solution for 5^ can be constructed by considering the terms with a decreasing 
number of Laplacians, and reads as 



= V2( ^^ a ^ +f(^ + ^ a )coth(Ct) 



5, £ + -f V'B j coth(Ct) (44) 



a 



+ \A a + ^V 2 Aj coth(Ct) 
The general solution of Eq. (J4"3~j) is 



& ~ o ^ 1 ^ ' 
where S °p is the general solution of the wave equation 



a 2 | ^ + 3ad 5 a p - V 2 5 °p = (46) 



^ ^ o 

satisfying S a/i „ = 0. 
o ^ 

Expanding the solutions of this equation in plane waves 

roo 

S Vx, t) = / D a Jk, t)exp(^k • x.)d 3 k , (47) 

JO 

the amplitude D a Jk., t) of the modes satisfies 

a 2 D a p + 3aoD a + (k • k)^ = (48) 

and is transverse, D a ^kp = 0. Changing to the variable x = sinh 2//3 (Ct) we get the following 
equation 

l(1 + I 3)|£ + ( 4l3 ^ + ^ , (49) 



where k 2 = f^c) ' k)> wr ^h the solution 



£> = x _3/ Vl + 4fc 2 a; + 4x 3 (50) 
x [P fc cos(a; fc (x)) + Q k sm(cj k (x))} 



and ( 

27 r x \[x' dx' 



/,« 2 < f x Vr dx' 
uj k (x) = \ k 6 H / ; . (51) 

V; V 64 io (x« + fcV + 1/4) vT+x 75 
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In summary, our solution for the linear perturbations of a k = pressure-free FRW universe 
with a cosmological constant has the form 

p = , 5p = , h 0a = — sinh- 2/3 (Ct)V 2 C a , (52) 

a 



= S a p + A af} coth(Ct) + 7] a pB (53) 

"^2V^ C ° th(Ct) + J 

8C 

+ — (C a ,fS + C^ a ) COth(Ct) , 

C M cosh(Ct) / 2 . 3V 2 5 r \ 

where A = A(x a ), B = B(x a ), C a = C a (x a ) are space functions, and the time-dependent 
amplitudes I and J are elliptic integrals, 

sinh 2/3 (Cr) . . 

J = -2- 1/3 3/(t) - 3sinh- 1/3 (Ct)cosh- 1 (Ct) . 

The term S a p represents gravitational waves and is a solution of the wave equation (J46|) . 

In obtaining this result we used the gauge freedom (jSJ) up to transformations with harmonic 
generator functions b and & satisfying 

V 2 6 = , VV = , = . 

By Eq.(|17|) the form of the solution is not preserved by arbitrary gauge transformations. 

4 The Sachs- Wolfe effect 

The Sachs- Wolfe effect is the contribution to the temperature variation ST of the cosmic back- 
ground radiat ion due to the gravita tional perturbations along the path of the photon. It can 



be computed ([Sachs h Wolfd. 119671 ) as follows, 



T 2 Jo \ Ot] Ot] 

where t]r and rjE denote the time of reception and emission, respectively, and w is the affine 
length along the null geodesic of propagation with tangent four-vector 

fH-i,o, < 57 > 

where e a is the three- vector of the photon direction normalized by e a e a = —1. We consider the 
contribution of the relatively increasing mode. Then /i 0| 8 = and the second term under the 
integral in Eq. (|5T)j) vanishes. The term S a p has the amplitude B a p. By using the relation 

()57|) in the identical decomposition 

y a,—, — dw = y a e a dw — y'dw (58) 
' dw 
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for an arbitrary function y = y(r), x a ), we get dipole anisotropy contributions with the amplitude 
B^e 13 and gr avitational redshift te rms. The time dependence of the trace part of h a p is different 
from that in ISachs fc Wolfel (|l967h thus the cancellation of the integrated terms does not occur 



here. Taken together the contributions to the temperature variation sum up to 

5T 3 riR-VE 



01 6 rVR-VE _ 



5 B, a e a A(t) - BA sw (t) 



2 7 / s a 2 C 2 



- VR 




+ 
- Ve 


(?) 



ISW 



at the respective events R of reception and E of emission, where the transfer functions of the 
fluctuations are defined 

A(t) = a(/coth(Ct))* = aosinh(Ct) 2/3 (/coth(Ct))* , (60) 
\sw(t) = yJ^c-2 S ™HCt)^A(t) (61) 



and 



(f ) = ^£lnh(a) 2 / 3 [5sinh(tt) 2 / 3 A(t)]^ 

1 /Q P R 

= -a C2 1/3 / BA ISW {t)dw (62) 

is the integrated Sachs- Wolfe (ISW) term with the transfer function 

= 2cosh 2 (Ct) + 3 J _ 2 1 / 3 
ISW 3sinh 2 (Ct) sinh 1/3 (Ct) cosh(Ct) ' 1 ' 

The terms containing B iCt e a have a dipole character (jSachs fc Wolfel . Il967t ). The relativistic 
redshift effect, represented by the second term in (jEHj) is displayed on Fig. 4. The amplitude 
Asw of the temperature variations decays exponentially. In the limit of a vanishing cosmo logical 
constant, the amplitude has the constant value 0.1. 



5 Physical interpretation 



The combined results of the WMAP and SDSS surveys (jTegmark et all l2004al ) are for the 
matter and dark energy densities of the Universe that fl m = 0.3 ± 0.04 and fl\ = 0.7 ± 0.04, 
respectively. For the Hubble parameter Ho = h ■ lOOkm/ sec/Mpc we use h = 0.7 ± 0.04. 
From these parameters we obtain the value of the cosmological constant A = 3HqQ\/c 2 = 



0.12 • 10 



-51 



m 



13.47 • 10 9 yr. The constant C can 
be determined by use of (fT3j) to be C = 2.85 • 10 _18 sec _1 . Taking the decoupling to occur at 
redshift z = 1100, this corresponds in the present model to td ec = 4.66 • 10 5 yrs. 



and the current age of the Universe tn = 

-18„„„-1 



Th e amplitude of the density contrast at the decoupling is estimated in Kolb fc Turner! 
(|1990h . to be 8p/ p ~ 10~ 3 . A more recent estimate based on the as measurements ( Tegmark et al 
l2004bl ) confirms this order of magnitude. The current value is determined by the amplitude of 
the relatively growing mode 




(/coth(Ct)) < 
(Icoth(Ct)) tc 




0.8 



(64) 
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This corresponds to an amplificat ion factor 800 . Our result is more than an order of magnitude 
higher then the data in Table 1 of lHeathl (|1977l ). lLahav et al.l ([19911 ) find that the growth factor 
depends only weakly on the value of the cosmological constant, but is much more sensitive to 
that value for higher redshifts. 

The amplitude of the temperature fluctuations 5T/T of the CMBR due to the Sachs- Wolfe 
effect on the relatively growing mode vanishes at Ct = 1.02. At the moment of reception, the 
amplitude is Asw — —0.033. Thus our conclusion is that the observed angular variation of the 
temperature fluctuations of the CMBR are attenuated by a factor of 3 in the presence of dark 
energy. The temporal va ri ation of the SW amplitude has been observed in detailed numerical 
studies ( Hu fc Sugivama ( 1995h ). These computations follow the variation of the multipole 
moments, which obscures the effect. 




Figure 4: The dependence of the transfer function Asw ° n the distance traveled by the photon from 
the event of last scattering. The decoupling can be taken to occur at Ctd ec = within the error bar. 
Reception is indicated at Ctn = 1.21. 



The weighing function of the ISW effect is displayed in Fig. 5. T he dominance of the 



late-t ime contribution is apparent, in agreement with earlier predictions ([Crittenden fc Turok . 
1996}). The correlation of the ISW powe r spec trum with the mass distribution from galaxy 
counts has been verified by iFosalba et al. (|2003h . 



6 Concluding remarks 

The results of our work are relevant in two major contexts: first, we find that the observed 
CMBR power spectrum is attenuated by a significant overall factor due to the presence of 
th e cosmological constant . Our finding is consistent with the results of the numerical studies 



m 



Hu fc Sugivama ( 1995| ). This phenomenon may help understanding structure formation by 



taking into account a higher level of the initial matter fluctuations. The second observation 
of interest that we make is the natural complexification of the perturbed fields in the normal 
picture of Legendre. This complex description may prove useful in a treatment of quantum 
fluctuations, a phenomenon that occurs also in the Fourier expansion of the wave solutions. 
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Figure 5: The weighing function of the ISW effect. 

7 Appendix: Legendre normal forms 

The integral ()29|) can be brought to the Legendre normal form (|Grobner &: Hofreiter . 1957 ) 



1 / i/o \ 1/2 3 + 3^0 
3 + 3 1/2 i) E F 



2V6 

2 1 / 3 (l-x)x 1 / 2 



6 



(65) 



(S 3 + 1)V2 

using the mirror symmetries of the incomplete elliptic integrals 



E = E 



( 




1/2 






2(as+l) 







and 

To get the Legendre form of 

J 



F = F 



( 


[■(a+aVaijaf 


1/2 






2(a:+l) 







3 /•* 



o w' 4 / 3 (l + w /2 ) 1/2 2 7o x' 3 / 2 (l + a;' 3 ) 1/2 



(66) 



(67) 



(68) 



we use the variable w = x 3 ^ 2 in the integral. The normal form is ( Grobner fc Hofreiterul957 ) 

j= - 3 2 i/2 ( 3 + 3 i/2 z) i/2 ( E _ 3 + 3l/2i F ) (69) 



(x 2 -x + l) l l 2 
3 (x + l) 1 / 2 ^/ 2 ' 

The following relation holds between the incomplete elliptic integrals: 

J + 



2V 3 ' 



a;i/ 2 (l + x 3 )i/ 2 



(70) 
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